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it tends to values larger than the free-pressure one. 


The convergence properties of the resummed thermal per¬ 
turbation series for the thermodynamic pressure are investi¬ 
gated by comparison with the exact results obtained in large- 
N cj>^ theory and possibilities for improvements are discussed. 
By going beyond conventional resummed perturbation the¬ 
ory, renormalization has to be carried out nonperturbatively 
yet consistently. This is exemplified in large-A (f>\ and in a 
special large-A t/ig model that mimics QED in the limit of 
large flavour number. 

I. INTRODUCTION 

A few years ago, the authors of Ref. Q have accom¬ 
plished the task of calculating all contributions to the 
thermodynamic pressure in QCD that are accessible to 
conventional resummed perturbation theory, that is up 
to contributions from the nonperturbative magnetic sec¬ 
tor which involve at least 5®. Dismayingly, the apparent 
convergence of the perturbation series up to order , 
rendered in Fig. la, proved to be too poor to warrant its 
usage at realistic couplings g « 2. In fact, the conver¬ 
gence is spoiled in particular by the contributions involv¬ 
ing resummation effects, as can be seen in Fig. la from 
the large jumps whenever a new contribution involving 
odd powers in g is included. 

Since this problem arises already before one has 
to face the nonperturbative nature of the (chromo)- 
magnetostatic sector, which contributes to 0 {g^T^), it 
may have to do with the constraint inherent in perturba¬ 
tion theory to truncate everything to finite-order polyno¬ 
mials in the coupling (up to possible logarithms). This 
is necessary in view of the needs of the renormalization 
process, but it certainly discards a large portion of the 
resummation effects which are in principle contained in 
low-order diagrams. 

The replacement of truncated perturbation series by 
Fade approximants, i.e. by perturbatively equivalent ra¬ 
tional functions, as proposed in Ref. Q can be considered 
as a very simple guess what the discarded terms may have 
been. At any rate, it gives a rough idea as to their im¬ 
portance. In Fig. lb the result for the pressure is trans¬ 
formed by appropriate Fade approximations, and it does 
have a big effect for larger coupling. The apparent con¬ 
vergence is somewhat improved, but at g ~ 2 the result is 
still inconclusive; actually, the final result including the 
5th order contributions looks even worse than before as 



(a) 



(b) 

FIG. 1. a) The perturbative results for the pressure of 
QCD(A/ = 3) up to order (full line). Short, medium, 
and long dashes give the results up to , g^, and g"^, respec¬ 
tively. b) The corresponding Fade approximants [0,2], [1,2], 
[2,2], and [3,2]. 

In this talk, I shall discuss the convergence properties 
of resummed thermal perturbation theory in simple solv¬ 
able but nontrivial cases. The starting point will be an 
exact nonperturbative formula for the full pressure which 
itself has the form of a one-loop integral. Its evaluation 
and renormalization will be carried out first in the large- 
N limit of scalar (()^-theory, which is simplified by the fact 
that the self-energy is momentum-independent, and then 
in a special large-(/)g model that mimics QED in the 
limit of large flavour number, which brings in some of the 
complications coming from momentum-dependent self¬ 
energies. Moreover, these toy models have some notewor¬ 
thy features which also make them interesting on their 
own. 
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II. A ONE-LOOP FORMULA FOR THE 
NONPERTURBATIVE PRESSURE 


III. NONPERTURBATIVE RENORMALIZATION 
IN THE EXAMPLE OF LARGE-A 01 THEORY 


In Ref. 1^ we found that at the cost of introducing a 
further integration with respect to the bare masses of a 
field theoretical model the full thermodynamic pressure 
can be expressed as a one-loop integral involving the full 
self-energy. This follows from the observation that 


where C is the Keldysh contour of real-time finite- 
temperature field theory HR and from the requirement 
that the pressure PiT) — P[0) 0 when all masses are 

sent to infinity.^ Thus 


PiT) - P{0) = -(t = 0^ 



<Pq 

(27r)" 


[£'t (9;too) 



cPq f 1-1- 2 n(( 7 °) 

(27r)" 1 _ jj^'2 _ Hriq^, q^; TOq) 


q^ — ttIq — n(q^; TOq) 


( 2 ) 


where denotes the real-time matrix propagator and 
n is obtained by diagonalization of the self-energy ma¬ 
trix. This can be easily generalized to fermions and in 
principle even to (gauge-fixed) gauge theories ||^. While 
it is true that the introduction of auxiliary masses gen¬ 
erally spoils both gauge and BRS symmetries, the com¬ 
binatorics of the loop expansion which is reshuffled by 
the mass integration in Eq. does not depend on these 
symmetries. So although their loss would be a high price 
in intermediate steps, the full integral depends only on 
the physical masses (which are zero for gauge bosons). 

One noteworthy aspect following from the essential 
one-loop nature of the above (exact) formula is that it 
appears to be manifestly infrared-finite in all dimensions 
n > 3, irrespective of whether or not IIt gives rise to 
screening (or even divergences). So this seems to be 
a good new starting point for some reorganization of 
thermal perturbation theory—which would have to see 
that when obtained in some approximation is not 
expanded out again from the denominator in Eq. i). 
Clearly, this would require a correspondingly reorganized 
renormalization procedure. So far, everything has been 
written down in terms of unrenormalized quantities only. 


'^For a concise review see the lectures of Peter Landshoff at 
this meeting 

^This approach has obvious similarities with the so-called 
exact renormalization-group approach of Ref. ^ 


The large-TV limit of a scalar 0(A) model with inter¬ 
action Lagrangian 


can be solved exactly It can also be viewed 

as a certain approximation to ordinary ^0‘^-theory, in 
which only those diagrams are kept which have a topol¬ 
ogy that corresponds to the leading term of a 1/A- 
expansion. This infinite set of diagrams is alternatively 
called Hartree-Eock, super-daisy, cactus, or foam 0 dia¬ 
grams. 

In this approximation, the Schwinger-Dyson equation 
for the self-energy does not involve vertex functions and 
is given by a one-loop equation. Its renormalization at 
T = 0 gives 


= rrig -I- AoM(to^) 


( 4 ) 


with 


M(to 2 ) = i 


(Pq i 

(27r)’^ q^ — m? -\- ie 


( 5 ) 


Similarly, coupling constant renormalization is given 
by 


A — Ag -f XqXM' ifn^) (6) 

when renormalizing at zero momentum. In four dimen¬ 
sions, this leads to the problem of triviality if one requires 
both Ag and A to be positive, for then A ^ 0 as n ^ 4. 
Keeping A > 0 is only possible by using a cut-off or by 
accepting Ag < 0. In the latter case one finds that the 
scattering amplitude has a tachyonic pole, but this oc¬ 
curs at s = Stach ~ ~ which is 

exponentially huge for reasonably small coupling. Either 
way we can accept this theory as an effective theory for 
momenta and temperatures <C v^tach- 

At finite temperature, the Schwinger-Dyson equation 
leads to a “gap” equation of the form 

-I- = TOq -I- AgM7’(TO^ -I- dm^) (7) 


with 


Mt(to^) - M(to^) = At(to^) 



Elimination of the unrenormalized parameters in 
favour of the renormalized ones yields 

Sm^inP, T) = A[M(m^, dm^) -|- Nxim^ + dm?)] (9) 

with 
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= M{m'^ + Sm^) — M{rn?‘) — (10) 


The function M is formally of order A^, i.e. of the same 
order as 3-loop contributions to the pressure, and it has 
been occasionally missed in the literature [|| . It exhibits 
a nontrivial interplay of the thermal mass correction 
6rn? with the zero-temperature UV-divergent quantities 
M{m?) and M^rn?) appearing in the counter-terms. 

Rewriting also the pressure in terms of renormalized 
quantities, we arrive at the remarkable formula 


P{T) - P(0) 





5rn?{m''^,T) 


( 11 ) 


where is a running coupling that is equal to A 

when m'^ = m?; Aq is kept hxed throughout. 

Fortunately, Eq. (|ll|) can be integrated, yielding 


P(T) - P(0) = + W) + 

1 X 

oo 

( 12 ) 


(see Fig. 2), which we have agreed to ignore. In order that 
•\/®tach be exponentially far away, we need A(T) ^ 10^, 
which also cuts out the case of no solution for the thermal 
mass and therefore for the pressure.^ 



FIG. 2. The relation between solutions for 5rr? and the 
renormalized conpling A; the heavy line marks where Sm 
would become equal to the (modulus of the) tachyon mass. 


IV. COMPARISON OF A NONPERTURBATIVE 
RESULT WITH PERTURBATION SERIES 


While the second term on the right-hand-side can be 
identihed as an interaction contribution, the subsequent 
terms again come from the thermal mass shift in zero- 
temperature integrals. 

While being UV-finite, the above equation becomes 
IR-divergent when m 0. This is not a problem specihc 
to hnite temperature, but comes from the breakdown of 
the on-shell renormalization scheme that we have used 
so far. In the limit m —> 0 one can switch to an off-shell 
(MS) scheme through 

which simplifies 

— 1 ^TTl ^ 

M (log - 1) (14) 

It is now only a matter of simple numerical integrations 
to calculate 6m‘^ and the resulting pressure. 

Somewhat surprisingly, there is not one but two solu¬ 
tions for Srri^ for small values of the coupling, and none 
for A > Acrit, whose value depends (in a renormalization- 
group invariant manner) on the renormalization scale fl. 
However, the higher values of dm? that correspond to one 
and the same A are close to the tachyonic scale -^/stach 


Besides its pedagogical value, the above nonperturba- 
tive result can be put to use to study the properties of 
(truncated) series expansions in X(U). The first few terms 
in the result for the pressure read^ 

P(T) - P(0) _ 7r2 A , A3/2 

T4 " 11 ^^ 576 ^ 

-(6 - 7 - log ^);^g432^2 

n A5/2 

-f(3 27 

- ((6 - 7 - log - 30 + ^) 

+0(W2) (15) 

Comparing truncated series expansions like the above 
with the full nonperturbative result, we can investigate 
the convergence properties of an expansion in powers 
(and logarithms) of A^/^. It turns out that these de¬ 
pend strongly on the ratio of the arbitrary renormaliza¬ 
tion scale p, to temperature. 

In Fig. 3 we juxtapose the exact and the perturbative 
result for the pressure [P(r) — P(0)] to its ideal-gas value 


^An explicit expression of the pressure in large-A cj)^ has 
been obtained previously in Ref. [^, though their formula 
does not satisfy the physically-important constraint that the 
pressure vanishes when the mass is infinite. 


^This is in accordance with the findings of Ref. 

®Up to and including order where there is no differ¬ 

ence between the foam-diagram subset and the full set of di¬ 
agrams, this agrees with the full N = 1 result obtained using 
resummed perturbation theory in Ref. 0- 
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7r^T'^/90, including successively up to 10 terms beyond 
the leading one. We choose various values of the renor¬ 
malization scale fl, but for ease of comparison in each 
case we plot against A evaluated for ft = T. 

This shows that when p, is very different from T, 
the convergence of the series deteriorates significantly. 
For p = lOOT (Fig. 3a), the truncated series develop 
oscillatory behaviour for larger values of the coupling, 
whereas for p = the perturbative results fail to 

improve with increasing order at roughly the same value 
of the coupling, albeit in a more peaceful manner. With 
p « 27rT, the perturbation series tolerates the largest 
coupling strength, but it is evident that the perturbative 
results cannot describe more than a few percent deviation 
from the free pressure value. 





^^^^^^_j 

1 2 3 4 5 6 7 8 

a1(T) 


FIG. 3. A comparison of the perturbative results for 
P = [P{T) — P{0)]/ '"^ as a function of A^^^(r) up to 12th 
order for different choices of the renormalization scale; a) 
p = lOOT, h) p = 2nT, c) p = 


In Fig. 4 the results of a Fade improvement is dis¬ 
played. This is seen to work surprisingly well except for 
those cases where the Fade approximant develops a pole 
beyond which the result appears to be off by a constant. 

In QCD, the Fade approximants did not work nearly 
as well (cf. Fig. lb). This might have to do with the ab¬ 
sence of log(A)-terms in the large-limit of our model. 
In QCD similar logarithmic terms are present and they 
are treated like constants when building the Fade ap¬ 
proximants. 



FIG. 4. As in fig. 3b, but with the perturbative power 
series replaced by Fade approximants [0,2], [1,2], [2,2], 
[2, 3],..., [4, 4]. They quickly converge to the exact result 
with the exception of [3, 3], which has a pole. 

Let us finally try out a completely different expansion. 
If we a posteriori expand everything in powers of 6 mjT 
rather than A, we find that the convergence properties 
are excellent so that substantial deviations from the free 
pressure results can be covered with only a few terms kept 
(Fig. 5). The reason for this is that the high-temperature 
expansion has a convergence radius of 6 m jT = 27r, a 
value which, at least in our model, is always far away for 
all A Acrit- 



FIG. 5. The pressure with solutions of the thermal mass 
equation truncated at leading order in 5m/T (top curve) and 
next-to-leading order (bottom curve). All higher approxima¬ 
tions are virtually indistinguishable from the exact solution, 
and they collectively form the middle line, which broadens 
only at the highest values of A. 
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This, perhaps merely tantalizingly, indicates that a re¬ 
organization of thermal perturbation theory as a series 
in 5m/T rather than g should lead to dramatic improve¬ 
ments. In this connection see also the proposal of a 
“screened perturbation theory” in Ref. |l^ and the con¬ 
tribution of S. Leupold at this conference |13|. 


V. NONPERTURBATIVE RENORMALIZATION 
IN THE EXAMPLE OF LARGE-A (f>l THEORY 


Large-iV 0^-theory was exceptionally simple in that the 
self-energy is momentum independent. In order to learn 
how to deal with the complications from a momentum- 
dependent self-energy and the need of wave function 
renormalization we now consider a scalar (^g theory, in 
which a scalar isovector A is coupled to N —>■ oo scalar 
isodoublets according to Q 

= r = l,...,N (16) 

so that this model mimics large-A/ QED, but without 
the complications of spinor and vector boson helds. The 
masses of the fields A and ip are denoted by m and M, 
respectively, and we require m < 2M in order to have 
stable photons. 

The same trick as in Sect. 2 can be used to derive a 
one-loop formula for the full pressure, but it suffices to 
differentiate and integrate with respect to the “photon” 
mass TO, for to —*■ oo switches off all interactions. Thus 


PiT) = P--(T.MS - |l„. Qra- j 

1 + 2n(g°) 1 

~ e- m'i - gliriq^) 

where 7r7’((7°,q^) and 7r(g^) are the finite-temperature 
and zero-temperature photon self-energies (with two 
powers of the coupling factored out), respectively. 

In the limit of A ^ oo, the interaction pressure is of 
order A°, and to this order there are no corrections to 
the internal “electron” lines in the one-loop photon self¬ 
energy, because the electron self-energy is of order A“^. 
So the latter can be ignored except in P|?^®®(T,Mg) ^ 
A, where the renormalization of Mg contributes to the 
interaction pressure, and is in fact essential to make the 
latter finite. 

The UV divergences in 7r(g^) have to be eliminated by 
mass and wave-function renormalization. Choosing on- 
shell renormalization we have 



to '2 


= rn'o -bgo7r(TO'^) 

= Z2{m'^)gl 


(18) 

(19) 


Z 2 {m'^) 

^T(9°,q^) 


- 1 _ dm'^ 

1 — (jfg7r'(TO'^) dm^ 

= Tr{q^) — — {q^ — m''^)TT' ( 21 ) 

= 7r'r(q°, q^) — 7r(TO'^) — {q^ — m''^)TT’{m'^) (22) 


In Eq. ©) the mass integration can be carried out, 
which is trivial before renormalization and only slightly 
less so after, with the result 


p(y)mt ^ 


g=0 


= -Im 

2 . 


d^q 

(27r)6 


2 n(g°)log 


g^7rT(g°,q^)+m^-g^ 
TO^ — 


g^7rT(g°,q^) -bTO^ - 
g'^Ttlq^) +m^ — q^ 

g^^T(g°,q^) ' 

g^7r(g2) -|- to^ — g^ 


(23) 


The last term in the square brackets comes from 
P^^^{T,Mq) — and is necessary to remove 

the UV divergence contained in the second term. There 
is a subtlety involved here in that ir = ttt — ^ except 
for spacelike momenta; for those this equality holds for 
the real parts only, while the imaginary parts differ by 
contributions that vanish exponentially in the ultraviolet. 
This again provides an example of a nontrivial effect of 
zero-temperature renormalization in thermal quantities. 

Beyond the large-A limit, the renormalization of 
Eq. becomes much more involved and requires 

a careful treatment of overlapping divergences, see 
Ref. 0. 

In the large-A limit, we had to consider only one- 
loop (but this time non-local!) contributions to the self¬ 
energies—still, Eq. (0) is clearly a nonperturbative re¬ 
sult involving arbitrarily high powers of g. 

It is in fact a result which could not have been obtained 
by the methods of hard-thermal-loop resummation. The 
hard-thermal-loop approximation to ttt turns out to be 
identical in form to the longitudinal component of 
in four-dimensional gauge theories, but with a reversed 
over-all sign: 

(9 > q ) = 

rp 2 

At g° this leads to a screening mass squared with a wrong 
sign, giving rise to a spacelike pole in the hard-thermal- 
loop resummed propagator. This has been noted before 
in (()g-theories , and it is similar to what occurs in the 
hard-thermal-loop graviton propagator when evaluated 
(inconsistently) on flat space [^, where this is inter¬ 
preted as the Jeans instability of a gravitating medium. 

In our scalar case, it is a reflection of the unbounded¬ 
ness of the potential from below. Eor small temperatures, 
the mass term for the photon stabilizes the system, but 
thermal mass corrections lead to a diminished total mass 
up to a point T = T^rit where the screening mass van¬ 
ishes. Beyond this point, an instability develops, because 
thermal fluctuations become able to surmount the bar¬ 
rier provided by the zero-temperature mass term in the 
potential, and this is reflected by the appearance of a 
“tachyonic” screening mass (Eig. 6). 


1 - 


1 90 , 90 + |q| 

2|q| 90 - |q| 


(24) 
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Correspondingly, the pressure (p^ ) is well-defined and 
real only for T < Tcrit- A hard-thermal-loop approxi¬ 
mation fails because it approaches the problem from the 
wrong “side”. 






FIG. 6. The location of poles in the “photon” propagator 
for g = 10 at T = m = mi = m 2 (a), which is below the criti¬ 
cal temperature, and at T = 1.5m (b), which is above it. The 
two dashed lines mark the light cone and the zero-temperature 
mass hyperboloid go = q^-fmi. The line between them is the 
location of the poles in the thermal propagator. In (b) the 
second full line below the light cone marks zeros of the real 
part of the inverse propagator, which correspond to poles of 
the propagator only at its intersections with qo = 0 because 
of large Landau damping for go > 0 and g^ < 0. 

On the other hand, the full nonperturbative N —>■ 00 
result (|^) can be evaluated numerically (if tediously) by 
a couple of nested numerical integrations {ttt cannot be 
expressed in terms of elementary functions, but has a 
rather involved analytic structure). 

In Fig. 7, the results of such a numerical evaluation is 
given for three values of the coupling g and m = M, and 
it is compared with the strictly perturbative g^-part of 
the potential (without hard-thermal-loop resummation). 

The nonpert urbat ive pressure exists up to the point 
T = Tcrit ~ •\/247rm/g. Presumably there is a finite 
but complex analytic continuation beyond this point, but 
our derivation is in terms of a manifestly real quantity. 


which becomes infrared singular for T > Tcrit. Right at 
the critical temperature, however, where (|^ ) is still well- 
defined, we are as close as possible to the behaviour of 
four-dimensional gauge theories in that there is a van¬ 
ishing screening mass, as is the case for magnetostatic 
modes in QED. 


P(TY^I(Cg'^T'^] 



FIG. 7. Numerical evaluation of the interaction pressure 
( p^ for three values of g in comparison with the leading per¬ 
turbative contribution. 


VI. CONCLUSION 

At least as concerns the thermodynamic pressure of 
hot QCD, conventional (hard-thermal-loop) resummed 
perturbation theory, even after including everything that 
can be treated perturbatively, falls fails to give reliable 
results for moderately large coupling. The perturbative 
series shows a surprisingly poor apparent convergence, 
far worse than usually the case for the range of coupling 
considered, with the consequence that only a few percent 
of deviation from the ideal-gas result can be described. 

It could well be that this bad behaviour comes from a 
certain incompleteness of the resummation of hard ther¬ 
mal loops, because finally everything has to be broken 
down to a truncated series, that is, a polynomial in g 
and log(g). 

Fade approximants are the simplest possibility to re¬ 
place the latter by functions that are more likely to fit 
the actual behaviour at larger values of g, but in QCD 
they only lead to a marginal improvement, if any. 

In order to help prepare the (uncertain) way to alter¬ 
native resummation schemes, we have considered sim¬ 
ple scalar theories. Going beyond a strictly perturba¬ 
tive scheme in the coupling involves a correspondingly 
nonperturbative renormalization. In the simple models 
that we have considered we have demonstrated the im¬ 
portance of a consistent renormalization scheme. 

In large-A^ (j)^ theories we have been able to compare 
thoroughly the convergence behaviour of an expansion 
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in the coupling with the full nonperturbative result. We 
have seen the importance of the choice of renormalization 
scale, but have found that the convergence of these ex¬ 
pansions is limited to small coupling and that the range 
of the latter increases rather slowly with the order of 
perturbation theory. In contrast to QCD, Fade approxi- 
mants worked surprisingly well, but this may be specific 
to the model considered, for which the pressure fortu¬ 
itously did not involve logarithms in the coupling con¬ 
stant. An exceedingly good approximation was obtained 
by an expansion in terms of 5m/T in place of g, but it is 
of course totally unclear how to implement such a scheme 
in more complicated theories. 
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